Exactly solvable spin ladder model with degenerate ferromagnetic and singlet states 
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I. INTRODUCTION 



D.V.Dmitriev ^ , V.Ya.Krivnov ^ and A.A.Ovchinnikov ^ 

Max-Planck-InsUtut fur Physik Komplexer Systeme, Nothmtzer Str. 38, 01187 Dresden, Germany 
and Joint Institute of Chemical Physics of RAS, Kosygin str. 4, 111911, Moscow, Russia. 

We study the spin ladder model with interactions between spins on neighboring rungs. The 
^\ ' model Hamiltonian with the exact singlet ground state degenerated with ferromagnetic state is 

0^ . obtained. The singlet ground state wave function has a special recurrent form and depends on two 

0^ ' parameters. Spin correlations in the singlet ground state show double-spiral structure with period 

of spirals equals to the system size. For special values of parameters they have exponential decay. 
^ , The spectrum of the model is gapless and there are asymptotically degenerated excited states for 

Q ' special values of parameters in the thermodynamic limit. 

(N 

There has been growing interest lately in quantum spin systems with frustrated interactions . Of special impor- 
■ tance are models for which it is possible to construct an exact ground state. The first example of such a model is 
the well-known Majumdar- Ghosh model ||^ which is the s = 1/2 chain with antiferromagnetic interactions Ji and J2 
of nearest neighbor and next-nearest neighbor spins, where J2 = Ji/2. Afterwards, a large class of ID models with 
exact ground state has been found and studied ]lO| . 

A considerable progress has been achieved in a construction of such models by using so-called matrix-product (MP) 
form of the ground state wave function |]ri|, f2|l. One of the example of the MP state is the ground state of the 
I Affleck-Kennedy-Lieb-Tasaki (AKLT) modelj^], which is the generic model of the Haldane phase of S = 1 chain. 
'"O ■ However, the MP wave functions have typically short-ranged correlations and, as a rule, describe systems with the 
^ ' gapped spectrum. 

In we have studied the model with exact singlet ground state degenerate with ferromagnetic state. This 

model has two different nearest neighbor and next-nearest neighbor interactions depending on the parameter v and 
is described by the Hamiltonian 
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with periodic boundary conditions and even N — 2M. 

In fact, this Hamiltonian describe the line of transition points from the ferromagnetic to the singlet state of the 
model 
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O ■ The ground state of is ferromagnetic (singlet) at (5 < ((5 > 0), where S = J13 + 2(1-%^) ■ When 6 = 0, the 
model ^ reduces to the Hamiltonian (^. 

It has been proved in that singlet ground state of (0) has zero energy at j/ > as well as the ferromagnetic 

^ , state, while all other states have the positive energies. It has also been shown that singlet ground state has double- 
spiral ordering (excluding some special values of parameter v). 

In we have also considered the second model which is equivalent to special case of the spin-^ ladder. This 
model depends on one parameter, has the non-degenerate singlet ground state and its ground state properties are 
similar to that of AKLT model. 

The exact singlet ground state wave function of these models has a special recurrent form: 

^q{M) = Po*M, (3) 

'^M = {sf + l^ist + l^2S^-- + 1^2sjf){st + Visj... + V2S^)-- 

x(4n-l+'^l4»-- + ^24)--(4-l+i'l4) lii ■•■ i) (4) 
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FIG. 1. The two- leg spin ladder. 

where sf is the s = ^ raising operator. Eq.(^ contains AI — operator multipUers and the vacuum state \ll ... J,) 
is the state with all spins pointing down. The function 'I'm is the eigenfunction of Sz with 5"^ = but it is not the 
eigenfunction of S^. Pq is a projector onto the singlet state. 

The models in [ p^ correspond to two particular cases of this wave function: vi — 1/2 — for the first model (|^), 
and i/i + + 1 = for the second model. 

In this paper we consider the general case of the wave function (|^) which allows us to construct new Hamiltonian 
depending on two parameters vi and 1^2 and having the degenerate singlet and ferromagnetic states. The singlet 
ground state of this Hamiltonian has double-spiral structure with period of spirals equals to the system size, but for 
special values of the parameters (which include the second model in |]l3|) the singlet ground state correlations have 
antiferromagnetic character with an exponential decay and the singlet wave function in these cases reduces to the MP 
form. The interesting feature of the model is the existence of the excited states which are asymptotically degenerated 
with the ground state for special values of parameters in the thermodynamic limit. 

The paper is organized as follows. In the Section 2 we will consider the construction of the Hamiltonian with the 
exact singlet ground state and will calculate spin correlation functions. In the Section 3 the Hamiltonian at special 
values of the parameters is considered. In the Section 4 the numerical results for the energy spectrum of the model 
are presented and Section 5 gives a brief summary. 

II. THE MODEL 

Now we will construct the Hamiltonian for which ^>a{M) is the exact ground state wave function. This Hamiltonian 
describes two-leg s = ^ ladder with periodic boundary conditions (Fig.l) and can be represented in a form 



(5) 



where describes the interaction between neighboring rungs. The spin space of two neighboring rungs consists 

of six multiplets: two singlet, three triplet and one quintet. At the same time, one can check that for open chain the 
wave function m contains only three of the six multiplets of each pair of neighboring rungs: one singlet, one triplet 
and one quintet. The specific form of the singlet and triplet components present in the wave function (^) depends 
on parameters vi and 1/2- The Hamiltonian /i^^i+i can be written as the sum of the projectors onto the three missing 
multiplets with arbitrary positive coefficients Ai, A2, A3: 
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where P^'^^^ is the projector onto the missing multiplets in the corresponding cell Hamiltonian. Actually, the Hamil- 
tonian contains also offdiagonal projectors between two missing triplets, but we use this freedom in advance to 
make exchange integrals on two legs of ladder and on each rung respectively equal, that is J12 — J34 and J13 = J24- 
The wave function (^) is an exact wave function of the ground state of the Hamiltonian with zero energy, 
because 
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i = 1,...M- 1 
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and Ai, A2, A3 are the excitation energies of the corresponding multiplets. 

So, '9 i\j is the exact ground state wave function with zero energy for the total Hamiltonian of an open ladder 
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Since the function m contains components with all possible values of total spin 5* (0 < < M), then the ground 
state of open ladder is multiple degenerate. But it can be proved by the same way as it was made in |l3yi4|] that for a 
cyclic ladder (||) only singlet and ferromagnetic components of 5* m have zero energy. Therefore, for a cyclic ladder (|^) 
'^q{M) is a singlet ground state wave function degenerated with ferromagnetic state. Besides, the following general 
statements are valid for the model (||) 

1) . The ground states of open ladder described by (||) in the sector with fixed total spin S are non-degenerate and 
their energies are zero. 

2) . For cyclic ladder the ground state in the S — Q sector is non-degenerate. The ground state energies for 
Q < S < M are positive. 

3) . The singlet ground state wave functions for open and cyclic ladders coincide with each other. 

Since the specific form of the existing and missing multiplets in the wave function (^) on each two nearest neighbor 
spin pairs depends on the parameters vi and 1/2, the projectors in (^) also depend on vi and 1^2- Each projector can 
be written in the form 

-Pfc'^ = •^12'' (^1 ■ S2 + S3 • S4) + j}3''(Si • S3 + S2 • S4) + j}4''Si • S4 + J2'^^2 ■ S3 

+ J^^\Si ■ S2)(S3 • S4) + ^Si • S3)(S2 ■ S4) + (Si • S4)(S2 • S3) + C« 

and this representation is unique for a fixed value of the parameters vi and V2- 

Substituting the above expressions for the projectors into Eq. we obtain the general form of the Hamiltonians 
hi^i^i. Inasmuch as the Hamiltonians ft-i.i+i have exactly the same form for any j, it suffices here to give the expression 
for hi^2'- 

hl,2 = Jl2{Ai2 + Asi) + Ji3(Ai3 + A24) + J14A14 + J23A23 

+J1A12A34 + J2A13A24 + J3A14A23 (10) 



where 



Aij — Si • Sj 



and all exchange integrals depend on the model parameters and the spectrum of excited states Ji 
Ji{vi, 1^2, Ai, A2, A3) as follows: 

A2 , X3 {1^1 ^ If - J^l , _ A2 A3 (J^i - 1)2 - 

■^12 — TT + TTT -, , 9 ■Jl3 — ' 
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z 

{Vl - lf{v2Vl +V2-I- vl) 

Z 

{Vl - lf[2vi - V2V\ - V2) 



z 



where 



^=J(^i-l)^ + J(^^i + 1)^(2^^2-^1-1)' 

(one should keep in mind that only positive A^ can be substituted to these expressions). 

In general, the Hamiltonian /ii,i+i contains all the terms presented in (|l^), but we can simplify it by setting, for 
example, J2 = J3 = and solving equations (O) for Ai, A2, A3. All Ai turn out to be positive in this case for any vx 
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and V2 except two lines: vi = \ and V2 = J^i + 1, where ground state is multiple degenerated. The Hamiltonian ft-i.i+i 
in this case takes the form 

hia = JuiAu + A34) + Ji3(^i3 + ^24) + JiiAu + J23A23 + JiAi2^34 (12) 
J12 - 2 -^13 2 

J/l(l - V2){'^1 - 1^2) 



Jii^vi-V2 Jta = - V2) Ji = 4- 



The calculation of the norm of (|^) and the singlet ground state correlation functions can be performed in complete 
analogy to the corresponding calculations for the case vi = 1^2 So, a norm of ^o(-^^) can be written in a form 



1 



Gm = (*o(M)*o(M)) = 2 y '^M{y)dy, (13) 



where $M(y) is expanded over Legendre polynomials Pn{y) 



M 



'i>M(2/) = ^c„(M)P„(2/) (14) 
The coefficients c„(^) are defined by the recurrent equation 

Cn[l + 1) = r ^ Cn-l{l) + Cn[l) 

In — 1 Z 1 

, n + 1 {V2(n + i)-v^-\f 

+2;rT3 2 (^^^ 

with initial condition co(0) ~ 1 and Cn(0 = at n > Z. 

The appropriate calculations result in the expression for spin correlation functions at ^ 00 

1 /47r/\ 

(SiS2, + i) = (S2S2,+2) = ^ COS — j (16) 

(SiS2,+2) = Jcos(^^ + A^) (17) 

These equations mean that the spiral on each leg with pitch angle ^ is formed and the shift angle between spirals 
on the upper and the lower legs is Ay; = ^ . So, there is just one full rotation of the spin over the length of the 
ladder, independent of the size of the system and for fixed I « N sX N ^ <x> two spins on the ladder are parallel. 

III. SPECIAL CASES 



There are special values of the parameters v\ and V2 for which Eqs.(16) and (|17|) are not valid. These values can be 
determined from Eq.(|l5|), when the coefficients of c„_i(Z) or c„+i(/) equal to zero. These conditions lead to following 
equations 

i^2(n-l) + !^i + l = (18) 



z/2(n + 1) + £^1 + 1 = 



(19) 



The special lines ( [1^ ) an d (|l9| ) on the (i^i, V2) plane arc shown in Fig. 2. At vx = V2 Eqs.([T^) and (19) give special 
points of the model (|ip ||l|,|14[. 

We note that there is a symmetry transformation 
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FIG. 2. {fi, V2) plane. Thick lines are special cases of the model. Each value of n corresponds to a pair of equivalent lines. 
Due to the symmetry (21) it is sufficient to consider only the shaded area of {v\, V2) plane. 



V2 
V2 - Vi 



(20) 



which maps lines ( p^ a nd ( |l9| ) onto each other and, therefore, there is pair of equivalent hnes for each value of n. 
The transformation~pO|) does not change the wave function (^) as well as the Hamiltonian (|l2|) apart from a change 
in energy scale by a factor (v^ — v-yf' . So, it is sufficient to consider the region on {1^1,^2) plane restricted by the 
inequality |i'2 — I'll < 1 (see Fig. 2). 

For vi and 1^2 defined by Eqs.(^8|) or ( p^ ^Miv) contains only n terms and the wave function ^'a/ contains only n 
multiplets rather than (M -I- 1) as it does in generic case. It can be shown jl^ that wave function ^'o(Af) for these 
special cases can be written in MP form: 

*o(M) = Tr {Di,2 D3A . . . Dn-i.n) , (21) 

where D ^ T + uS is the n x n matrix describing states of corresponding spin pair. Singlet state matrix is 

S^I \s) (22) 

where / is identity matrix and \s) is the singlet state. Triplet state matrix T is expressed by Clebsch-Gordan coefficients 

a 



nil ,m2 



((1, mi) (j, 7712) I (j, mi + 7712)) as follows: 



T = 



/ Co,, |0) Ci,,„i|l) 

' 1-1) Co.,-i|0) . 

... 

. . Ci._,|l) 

V 1-1) Co,-, |0) / 



(23) 



where j — ^^^^ and \a) is the triplet state with = a. The parameter u is defined by expression 
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V2 {n - 1) 

Exact calculation of the correlators in these cases using standard transfer matrix technique results in 



(24) 
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FIG. 3. Stripe spin structure on the ladder model. 
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) =-(l + 2z) 
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(25) 



where we use notations 



n - 1 



0^1 = 1 + 2z + U , UJ2 = LOl — Az 



In the particular case of zero singlet weight, u = 0, when spins on each rung form a local triplet, correlation 
functions ( psf ) coincide with those obtained in [Tl|Jl^ . 

According to Eqs.(|25|) the singlet ground state has coUinear or stripe spin structure, i.e. spin-spin correlations are 
ferromagnetic along legs and antiferromagnetic between them (Fig. 3). 

These correlations have an exponential decay for finite value n and the correlation length Vc is 



21n"^ 



UJ2 



(26) 



For n — > 00 there is a magnetic order m: 

(SiSi+2/) = — (SiSi+1+2/) 



I 



When It = ±1 the magnetic order is equal to the classical value 1/2. 

It is interesting to note that Tc diverges for finite n when u ^ 00 {{vi,V2) (~liO)) but the prefactors in ( |25| ) 
vanish in this case. As it is shown on Fig. 2 this is the point where all special lines are intersected and the wave 
fimction 5*0 (M) is a simple product of the singlet pairs 



4) lii - i) 



(27) 



We note that the wave function (^) show double-spiral ordering for all values of vi and 1/2 excluding the special 
lines. The crossover between spiral and strip states occurs in the exponentially small (at N 00) vicinity of the 
special lines. 

Now we should make the following remark. There is one particular case vi + 1^2 + ^ = {n = 2), which was 
considered in jl^ , when on each two neighboring rungs the wave function ^/m contains only one singlet and one 
triplet and does not contain quintet. Therefore, in this case the cell Hamiltonian can be written in the form jl^] 



(28) 



k=l 



where P, 



is a projector onto the quintet state. If Ai, A2, A3, A4 > 0, wave function 5'o(M) is non-degenerate 
singlet ground state wave function for the Hamiltonian (p8|). In this case all four-spin interactions can be excluded 
by setting Ji = J2 = J3 = and then we arrive at the Hamiltonian 



J12 



(S2i-1 



1 



) + (S2'i+1 • S2i 



+2 



+ 2Jl3 



(29) 



-f-2 Ji4(S2i-l • 82i+2 ^ ^) + 2J23(82i • S2i+1 — -) 
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and all exchange integrals Jij depend on one model parameter. The explicit form of Jy was found in It was 
shown that for this model \E'o(M) is non-degenerate singlet ground state wave function with exponentially decaying 
spin correlations and there is an energy gap. 

In other special cases we can also construct Hamiltonians for which ^'o(M) is non-degenerate singlet ground state 
wave function. But we have to introduce more distant interactions. For example, such a model for n = 3 would 
contain interactions between next-nearest neighbor spin pairs 

M 

H = His+i,i+2 (30) 

For the point [vi — \, 1^2 — when spins on each rung form a local triplet, one of the possible Hamiltonians 
can be written as 

-ffi,2,3 = A12 + A23 + 38Ai3 + 6^23 - {Al^ + Al^){Ai^ + J) 

-{Ar, + \)(Al^+Al^)^b2 (31) 

where 

Aij — Li • Lj — 1 

and Li = S2i_i -|- S2i is the S = \ operator. 

The ground state wave function for the Hamiltonian (|T]) has the matrix-product form ( pT|) with 

D = 




IV. SPECTRUM OF THE MODEL 

Generally, the excitation spectrum of the model (||, ^2|) can not be calculated exactly. It is clear that this spectrum 
is gapless because, for example, the one-magnon energy is ^ N^^ at N ^ 00. Moreover, the lowest singlet excitation 
is gapless as well. For the model (||, |l^) lying on the special lines this fact can be established from the following 
consideration. For the simplicity we consider the case 1^1 = 1^2 = = ^- We choose the variational wave function of 
the excited singlet state at i^i = = as 

^sii^, S, M) = \ [c{5, M) vE-ol^., M) + 5, M)] , 

where 4*0 (i^, M) and "^0(1^ + M) arc the normalized singlet ground state wave functions of (^) with zero energies at 
vi — V2 — V and vi — V2 — v -'r 5 (the point vi — V2 — i' + 5 on (i^i, V2) plane does not belong to the special line). 
The functions s{v, S, M) and '^o{i', M) are orthogonal and c((5, M) is the overlap of ^'o(i^, M) and "^{^{v + 5, M) 

c{5, M) = (*o(i^, M) \'fo{i^ + S, M)) 

It can be shown that 

^'('^'^^)- l + ^2(M!)2,o(M) «^ (32) 

Eq.(^2|) follows from the fact that (5'o(i^, M) \ '^o{v, M)) and c(i5, M) is represented by the sum of n terms in ([l^ ) 
while {^i^{v + 5,M) \ ^q{v + 5,M)) contains M terms. 

The variational energy calculated with respect to s{v, (5, M) is 
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1 



1 -c2(J,Af) 



{^o{'^ + S,M) I H{v,M) I *o('^ + '^,M)) 



1 + (52(M!)2 

(5(Af!)2 

And minimization of Es over the function 5{M) leads to 

J ^ g-Mln Af+0(J\/) ^ ^ g-A/lnAf+0(Af) 



This consideration can be easily extended to the points vi ^ V2 on the special lines. But it is not valid for the 
parameters vi and V2 which are out of special lines. We performed the numerical diagonalization of finite ladders for 
various parameters vi and V2 ■ The energies of lowest singlet and triplet states of (^) are shown in Fig. 4 as a function of 
iV for parameters corresponding to different types of the ground state. Figure 4 shows that the exponential degeneracy 
possibly takes place for all parameters vi and ^2, but we can not confirm it strictly. 

Thus, on the special lines the ground state of considered model is asymptotically degenerated at the thermodynamic 
limit. It is not clear if the degeneracy is exponentially large or not. 

So far, we have considered the models with degenerate singlet and ferromagnetic states. Now we discuss the phase 
diagram of the zigzag chain model given by (^) . The line of transition points from the ferromagnetic to singlet state 
is described by the one parametric Hamiltonian (|l|) (or by ( |l^ aX vi ~ V2) ■ 

The exact ground state in the singlet phase (Fig. 5) is generally unknown. But it is interesting to note that the 
groimd state on the line J13 = —1/2 is the product of singlets on ladder diagonals (2,3), (4,5), ... as in the point 
J13 = — J23 = —1/2 on the transition line. The spectrum of (|^) on the transition line is gapless. There are some 
regions on the plane (J13, J23) which were studied by different approximation methods. 

At Ji3 = and < J23 < 1 the model (||) reduces to the alternating Heisenberg chain studied in [0. The lowest 
excitation is the triplet and there is the gap. At J23 = and J13 > the model (||) reduces to the spin ladder with 
antiferromagnetic interactions along legs and the ferromagnetic interactions on rungs. It is evident that there is a gap 
at Ji3 ^ 1 (in this case the model is equivalent to the spin S ~ \ Heisenberg chain). It was shown in |jl^ that the 
gap exists at J13 ^ 1. At J23 = — 1 and J13 ^ 1 the spectrum is gapless according to the results of [ p^ . 

We have calculated the first singlet and triplet excitation at J13 = —1/2 and 1/2 < J23 < 1 by the numerical 
diagonalization of the finite ladders. As it can be seen from Fig. 6 the gap is closed on the transition line. So, we 
expect that the phase diagram of the model (||) has the form shown in Fig. 5. 
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-1 -0.5 



FIG. 5. Phase diagram of the zigzag chain model (2). The thick solid line is the boundary between the ferromagnetic and 
singlet phases. Circles correspond to the special cases of the model. The thin solid line denotes the heuristic boundary between 
gapped and gapless phases. On the dotted line the ground state is a product of singlet pairs [2, 3][4, 5].... 




FIG. 6. Dependence of singlet-singlet and singlet-triplet energy gap on J23 along the line J13 = —1/2 (dotted line in Fig.5). 
Calculation was made for finite chain with N — 20. 
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V. SUMMARY 



We have constructed the spin ladder model with ferro- and antiferromagnetic interactions between spins on neigh- 
boring rungs. The model has exact singlet ground state degenerated with ferromagnetic state. The spin correlators 
in the singlet ground state show double-spiral ordering with period of spirals equals to the system size. However, for 
special values of the parameters spin correlators in the singlet state have exponential decay and in these cases the 
singlet ground state wave function can be represented in the MP form. The spectrum of the model is gapless and 
there is asymptotic degeneracy of the ground state for special values of the parameters at the thermodynamic limit. 

The singlet ground state wave function has the recurrent form (|^) and depends on two parameters. This function 
can be further generalized. For example, we can take ^o(-^) in a form 

*o(M) = ^0*M 
where ^'m is the product of alternating multipliers 

I'm = (S^ + ViS^ + V2st-){st, + V^isi + V^st, ...){st, + VlS+ + V2S^. ..)... \ii ... i) 

with the condition 

1 + Vl _ 1 + 1^3 

The Hamiltonian of the model for which '^^{M) is the singlet ground state has two rungs in the elementary cell. 
This Hamiltonian can be taken in the form containing the interactions between neighboring rungs only and without 
four spin terms. 
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